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Abstract. Hadrons carrying heavy quarks, i.e. charm or bottom, are important probes 
of the hot and dense medium created in relativistic heavy ion collisions. Heavy quark- 
antiquark pairs are mainly produced in initial hard scattering processes of partons. While 
some of the produced pairs form bound quarkonia, the vast majority hadronize into 
particles carrying open heavy flavor. Heavy quark production has been studied by the 
PHENIX experiment at RHIC via measurements of single leptons from semi-leptonic 
decays in both the electron channel at mid-rapidity and in the muon channel at forward 
rapidity. A large suppression and azimuthal anisotropy of single electrons have been ob¬ 
served in Au + Au collisions at 200 GeV. These results suggest a large energy loss and 
flow of heavy quarks in the hot, dense matter. The PHENIX experiment has also mea¬ 
sured J/ i/r production at 200 GeV in p + p, d + Au, Cu -l- Cu and Au + Au collisions, both 
at mid- and forward-rapidities, and additionally Cu -l- Au and U + U at forward-rapidities. 
In the most energetic collisions, more suppression is observed at forward rapidity than at 
central rapidity. This can be interpreted either as a sign of quark recombination, or as a 
hint of additional cold nuclear matter effects. The centrality dependence of nuclear mod¬ 
ification factor, Raa(Pt), for J/il/in U + U collisions at = 193 GeV shows a similar 
trend to the lighter systems, Au -l- Au and Cu -l- Cu, at similar energy 200 GeV. 


1 Introduction 


The measurement of inclusive hadron yields in central Auh-Au collisions at RHIC led to the discovery 
of the suppression of hadron production at large transverse momenta (pj) compared to p + p colli¬ 
sions 111. This is generally attributed to the energy loss of light partons in the dense nuclear matter 
created at RHIC. Heavy quarks, i.e. charm and beauty, are believed to be mostly created in initial hard 
scattering processes of partons El and thus are excellent probes of the hot and dense matter formed in 
nucleus-nucleus collisions at high energy. While some of the produced pairs form bound quarkonia, 
the vast majority hadronize into hadrons carrying open heavy flavor. They interact with the medium 
and are expected to be sensitive to its energy density through the mechanism of parton energy loss. 
Due to the large mass of heavy quarks the suppression of small angle gluon radiation should reduce 
their energy loss, and consequently any suppression of heavy-quark mesons like D and B mesons at 
high-pr is expected to be smaller than that observed for hadrons consisting of light quarks 0. 
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We quantify the medium effects on high-p^ production in nucleus-nucleus collisions, A H- A, with 
the nuclear modification factor which is defined as following: 


Raa(Pt) - 


1 

(Neoll) 


yield per A -i- A collision 
yield per p + p collision 


1 d^N* * /dydpx 

(N,o„) ^ d2N'’"7dydpT 


( 1 ) 


This factor reflects the deviation of measured distributions of nucleus-nucleus, A -H A, transverse mo¬ 
mentum, at given impact parameter b, from measured distributions of an incoherent superposition 
of nucleon-nucleon (p + p) transverse momentum, scaled by the average number of expected binary 
collisions (Acoii )■ This normalization often is known as “binary collisions scaling". In the absence 
of any modifications due to the ‘embedding’ of elementary collisions in a nuclear collision, we ex¬ 
pect Raa = 1 at high-p.p. At low p^, where particle production follows a scaling with the number of 
participants, the above definition of Raa leads to Raa < 1 for < 2 GeV/c. 

However, the interpretation of heavy-ion collision data is complicated by the fact that heavy quark 
production is modified in a nuclear target by cold nuclear matter (CNM) effects. These include mod¬ 
ification of the parton density functions in a nucleus, parton energy loss in CMN, and (for quarkonia 
states), breakup due to collisions with nucleons in the target. To disentangle such effects, multi¬ 
ple measurements are needed, such as collisions of different ions and at different collision energies. 
RHIC has the unique ability to store and collide nuclei of different mass, i.e. Au -H Au, Cu -H Cu and 
U + U, and also collide asymmetric systems such asd + Au, He -H Au and Cu -H Au. By observing the 
modification in more elementary systems like in d + Au and He -H Au, one may hope to understand 
the initial modifications before the QGP is formed. 

Heavy quark production has been studied by the PHENIX experiment at RHIC via measurements 
of single leptons from semi-leptonic decays in both the electron channel at mid-rapidity and in the 
muon channel at forward-rapidity. In this paper I will summarize the latest PHENIX results concern¬ 
ing open and closed heavy flavor production as a function of beam energy and systems size. 


2 PHENIX Experiment 

The PHENIX detector comprises three separate spectrometers in three pseudorapidity ( 77 ) ranges. 
Two central arms at midrapidity cover | 77 |< 0.35 and have an azimuthal coverage {(p) of njl rad 
each, while muon arms at backward and forward rapidity cover -2.2 < 77 < -1.2 and 1.2 < 77 < 2.4, 
respectively, with full azimuthal coverage. In the central arms, heavy quark production was studied 
via measurements of single leptons (electrons) from semi-leptonic decays. Charged particle tracks 
are reconstructed using the drift chamber and pad chambers. Electron candidates are selected by 
matching charged tracks to hits in the Ring Imaging Cherenkov (RICH) counters and clusters in the 
Electromagnetic Calorimeter (EMCal). At forward and backward rapidity, the heavy quark production 
is measured via dimuon decays. Muons are identified by matching tracks measured in cathode-strip 
chambers, referred to as the muon tracker (MuTr), to hits in alternating planes of larocci tubes and 
steel absorbers, referred to as the muon identifier (MuID). Each muon arm is located behind a thick 
copper and iron absorber that is meant to stop most hadrons produced during the collisions, so that 
the detected muons must penetrate 8 to 11 hadronic interaction lengths of material in total. Beam 
interactions are selected with a minimum-bias (MB) trigger requiring at least one hit in each of the 
two beam-beam counters (BBCs) located at positive and negative pseudorapidity 3 <| 77 |< 3.9. 


3 Open heavy flavor - single electron/muon measurements 

Open heavy flavor production is measured in PHENIX through the measurement of inclusive electrons 
or muons Em. These analyses use a cocktail method to remove fake and real electrons/muons from 
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(a) e* and ifi : Au + Au vs. d + Au at Ws = 200 GeV (b) e* „ and ;r® : Cu + Cu at = 62.4 GeV 



Figure 1. The nuclear modification factors Rjaii(Pt;) and RAA(Pr) for th® ^hf ™d Pjjf '■ panels (a) and (b) for 
minimum bias d + Au, Au + Au and panel (c) for central collisions Cu + Cu. The two boxes on the right side of 
the plot represent the global uncertainties in the d + Au, Au + Au (Cu + Cu) values of Acoii (THB). 


the data sample, i.e. for the electron measurement, the electrons that come from either meson decay 
or photon conversions are measured and subtracted from the inclusive spectrum and the remainder is 
attributed to electrons coming from the semi-leptonic decay of D and B mesons. This remainder is 
also referred to as the non-photonic electron (NPE) component. PHENIX has measured spectra of the 
single electrons CHini and single muons MM from heavy flavor in p + p and Cu -H Cu collisions 
as well single electrons from heavy flavor in Au -H Au lfT^[T3]| and dn-Au IIT4HT51 collisions. 

The effects of CMN are expected to be present in the initial state of A -H A collisions; however, 
this CNM enhancement is convolved with the suppressing effects of hot nuclear matter. Eigure |l(a)| 
shows RdAuiPj) and Raa{Pt) for and the neutral pion for which only small CNM effects are 
observed between and for Pj <5 GeV/c. Above Pj ~5 GeV/c, where the CNM effects on both 
species, and are small, their Raa{Pt) values are consistent within uncertainties. However, in 
the range where CNM enhancement is large for and small on tt®, the corresponding efjp Raa{Pt) 
values are consistently above the values. This could suggest that the difference in the initial state 
cold nuclear matter effects due to the mass-dependent Cronin enhancement is reflected in the final 
state spectra of these particles in Au -i- Au collisions, although alternate explanations involving mass- 
dependent partonic energy loss in the hot medium are not ruled out. 
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Figure 2. The upper panel shows the nuclear mod¬ 
ification factor, as a function of centrality. The 
lines show theoretical calculations: the effect of cold 
nuclear matter (dotted, uppermost pair), direct pro¬ 
duction (dashed), and coalescence (dot-dashed). The 
solid lines represent the sum of all effects considered. 
The lower panel shows the ratio of forward- to mid¬ 
rapidity Raa Ref. m. 
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Figure 3. Nuclear modification factor, Raa as a 
function of centrality in Cu-l-Au collisions at forward 
(red) and backward (blue) rapidity. For comparison, 
the Au-)-Au (green) data are also shown, averaged 
over forward and backward rapidities Ref. ISD 


To study the interplay between initial-state effects and final-states effects for heavy flavor 
productions, PHENIX measured the Raa{Pt) in Au -H Au collisions at low beam energy, - 

62.4 GeV. These Raa{Pt) values for in Au + Au collisions at 62.4 GeV are compared to Raa{Pt) 
at the same energy as shown in Figure 1(b) At 62.4 GeV, the competition, if present, favors heavy- 
flavor enhancement over suppression. This is consistent with previous results with hadrons where the 
Cronin enhancement increases as the collision energy decreases 0. 

PHENIX measurements of heavy flavor muons ip^p) at forward rapidity (1.4 < y < 1.9) in central 
Cu-hCu collisions at - 62.4 GeV show a significant suppression 0. The magnitude of this 
suppression at forward-rapidity in Cu + Cu (shown in Fig 1(c) i is compared to the suppression of ef^p 
in central Auh-Au collisions at mid-rapidity at the same energy. As open heavy flavor is significantly 
more suppressed at forward-rapidity than at mid-rapidity in Cu -H Cu, additional nuclear effects, such 
as gluon shadowing at low-Bj6rken-.r or partonic energy loss in the nucleus, may be significant. The 
heavy flavor e* ^ and p^p are compared in Fig|l(c)|to a theoretical prediction that combines the effects 
of partonic energy loss, energy loss from fragmentation and dissociation, and includes nuclear matter 
effects such as shadowing and Cronin enhancement due to parton scattering in the nucleus im. While 
consistent within uncertainties, the model predicts more suppression for heavy flavor electrons than 
seen in the data. 


4 Heavy quarkonia - J/ij/ production versus system size 

Heavy quarkonia has long been proposed as a sensitive probe of the color screening length and de- 
confinement in the quark-gluon plasma ifTTll . It was originally suggested by Matsui and Satz ifTSl that 
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Figure 4. Nuclear modification factor, Raa(Pi}, measured in U + U (blue) collisions at =193 GeV as 
a function of centrality. For comparison, Au+Au (red) and Cu+Cu (green) data are also shown. The U + U 
Raa(Pt:) results are preliminary. 


a signature of the formation of the QGP would be that 7/i/r states that are initially bound may disas¬ 
sociate due to color screening in the high temperature QGP matter. The picture that was originally 
proposed is complicated by competing effects, which modify quarkonia production and survival in 
cold and hot nuclear matter. Heavy quarkonia measurements are made in PHENIX by either detecting 
opposite sign electrons at mid-rapidity (|y|<0.35) or by detecting opposite sign muons at forward- 
rapidity (1.2 <|y|< 2.2), reconstructing the invariant mass of the di-lepton pair, and subtracting the 
continuum background HU. The versatility of RHIC to provide heavy ion collisions at different cen¬ 
ter of mass energies and for different ion species has been further demonstrated during RUN-12 by 
providing first U + U and Cu -H Au collisions at =193 and 200 GeV, respectively. Both systems 
serve as an important test of the initial geometry. 

4.1 7/i/f in Au-H Au collisions 

The nuclear modification factor, RAAiPj), for Jlij/asa function of centrality (Npan) at mid-rapidity and 
forward-rapidity from Au -H Au collisions at aJJ^ - 200 GeV is shown in Fig. The data show that 
the suppression of J/if/ at forward-rapidity is stronger than at mid-rapidity, even though the energy 
density is expected to be slightly larger at mid-rapidity. This could be due to stronger cold nuclear 
matter effects at forward-rapidity, or possibly stronger coalescence of charm quarks at hadronization 
at mid-rapidity, or both. 

The comparison of the experimental data to the most recent theoretical calculations that incorpo¬ 
rate a variety of physics mechanisms including gluon saturation, gluon shadowing, initial-state parton 
energy loss, cold nuclear matter breakup, color screening, and charm recombination are presented on 
Fig. and described in details by PHENIX in Ref. ||T91 . 

4.2 J|4^ in Cu -H Au collisions 

Hot matter effects and CNM effects are present together in heavy ion collisions, and both are im¬ 
portant. In Au-hAu collisions at RHIC, for example, the addition of hot matter effects increases the 
suppression of the Jjxjj by a factor of roughly two over what would be expected if only CNM effects 
were present IT^EO). By colliding ions different sizes, like Cu -H Au, the cold nuclear matter effects 
is expected to be different at forward and backward rapidities. The comparison of d-i-Au, Auh-Au and 
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Cu+Au y/t/r modifications across rapidities may provide key insight on the balance of cold and hot 
nuclear matter effects, and whether they are factorizable. 

Figure 1^ shows nuclear modification factor, RAAiPj), for Jj^i as a function of Npan for Cu + Au 
collisions at - 200 GeV lISTl . The RAAiPj) for Au + Au collisions at the same collision energy 
and rapidity is shown in Fig. [^for comparison. The dependence of the Cu + Au nuclear modification 
on Npart at backward (Au-going) rapidity is similar to that for Au + Au collisions, while the Cu + Au 
Raa(Pt) at forward (Cu-going) rapidity is little lower. 

4.3 J|4^ \nU + U collisions 

The study of collisions of deformed nuclei like was initially proposed ca because they promise 
an additional gain of initial energy density relative to collisions of spherical nuclei. A significantly 
deformed initial geometry at very high energy density for specific orientations of U + U collisions 
is expected to have observable effects on elliptic fow, jet quenching, Jj^ suppression and other ob¬ 
servables that characterize the properties of the quark gluon plasma (QGP) |l22l|23]. The interesting 
configurations with high energy density can be selected in experiments by using a combined cut on 
the elliptic how value and the number of spectators for high multiplicity events El. 

The nuclear modification factor, Raa{Pt), for Jjxjj as a function of centrality, Npart, at forward- 
rapidity (1.2 <|y|< 2.2) from U + U collisions is shown in Fig. The Raa{Pt:) for Au -H Au and 
Cu -H Cu collisions at similar energy, ^[s^ - 200 GeV, and rapidity are also shown in Fig |4 for 
comparison. The centrality dependence ofU + U Raa{Pt) for J/if/ shows a similar trend to the lighter 
systems, Au -i- Au and Cu -H Cu. 


5 Summary 

PHENIX has shown first measurements at RHIC from two very unique systems: U + U and Cu -H Au 
at ^fs^ =193 and 200 GeV, respectively. Data from both systems were compared to Au + Au and 
Cu -t- Cu collisions at similar energy, - 200 GeV. The comparison aimed to shed more light on 
the role of the initial collision configuration and the final data observables. The U + U system holds 
the highest energy density so far measured at RHIC. We observed that the dependence of the Cu -H Au 
nuclear modification for J/if/ as function of centrality, Npart, at backward (Au-going) rapidity is similar 
to that for Au -i- Au collisions, while the Cu -i- Au Raa{Pt) at forward (Cu-going) rapidity is noticeably 
lower. Moreover, the centrality dependence of U + U RaaCPt) fo^^ -^1^ shows a similar trend to the 
lighter systems, Au -H Au and Cu -i- Cu at similar energy. These observation lead to the conclusion 
that the centrality dependence of Raa{Pt) for Jlijj at given energy is less sensitive to the system size. 
Therefore these new PHENIX results of Raa{Pt) for Jlijj 'mU + U and Cu + Au collisions are crucial 
for testing different theoretical approaches. 
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